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Abstract 


mspa.per presents an improved dynamic analysis for liquid annular seals with arhUmr,, 
profile based on a method , first proposed by Nelson and Nauven An j £ , ^ 

solution that incorporates a continuous interpolation of perturbed Quantities In /if*** ° ^ 
ferential dsrec„o n u Panted. The ongmaTmethcd uieTanZ^ZZ^slZeZZ- 
cum ferential gmdients, based on Fast Fourier Tnsnsforms (FFT). A simpler .Zeme baled 
hmhlZrr 3P r e ’, U , °T'‘ ‘° * com P uta ‘ wr ‘sll ! / more efficient with better convergence at 

varying seal profile in both axial and circumferential diZrZZ. V f ? arbltra r nl y 
elliptical seal with varying degrees o, aZZlZlm fsTZyZ A ZT&XL2 Tn 
actual operating clearances (6 axial planes with 68 cleamnces/plane) of an interstage seal 

prennflT M, “ n e ' gk PVe, ‘ ° Xmn P (SSME.ATD-HPOTP) ,s 


113 


NOMENCLATURE 


, 6 t 

<4 HI -4 Ml 

B n y B„, Bh 


cq 

Ci , c.. 

Cr,C } , 

C TT> Cyy 

Cflj' L 'yT 

c(z,J) 

F E 

frO' f,0 

f x :f u 

fx, fy 

ho, uo* I'o. Pq 

hi, u x , 1’i.pi 

h 

All. A yy 

I* Jt 

a iy » ^yar 

L 

M xt 

m 
POi 
Poe 
p$r 

R 

t 


xy i m yx 


li, u 


w 

\V 

Xq,Yq 

z. 3 

P 

P 

6 


0 


a? 


spatially dependent parts of first order solution 
coefficients of the variables of first order axial momentum 

equation _ , . 

coefficients of the variables of first order circumferential momentum 

equation 

nominal clearance (m) 

inlet and exit clearances (m) 

x, y axis clearances of elliptical seal (m) 

direct damping coefficients (N-s/m) 

cross coupled damping coefficients (N-s/m ) 

clearance function 

eccentricities along x and y axes (m) 
friction coefficients (Moody’s or Hirs ) 
x and y components of seal force (N) 
unbalance forces (N) 

variables of zeroth order (steady state) equations 

variables of first order (perturbed) equations 

film thickness (m) 

direct stiffness coefficients (N/m) 

cross coupled stiffness coefficients (N/m) 

length of the seal (m) 

direct mass coefficients (kg) 

cross coupled mass coefficients (kg) 

entrance pressure (Pa) 

exit pressure (Pa) 

pre-swirl ratio 

radius of the rotor (m) 

time (s) 

axial and tangential velocities (m/s) 
rotor surface velocity, u)R (m/s) 
preload 

axes of the elliptical whirl orbit 
axial and circumferential coordinates 

density (kg/m 3 ) 
dynamic viscosity (Pa-s) 
ellipticity, (c x - c y )/c x 
eccentricity ratios 
external load angle (rad) 
entrance loss coefficient 
angular frequency (rad/s) 
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INTRODUCTION 

sure^ Oxy^e^TWb^ump *"( A^D^HP OTTP ^T C ShU “ leM f m En ^ me ( SSM E) High Pres- 

investigated utihzing finite element models’^ thTent^pu™ ^ZTl ^ ^ 
designed with either straight or t^n^r^A c , P um P* Annular seals, initially 

distorted during the CsfoHM^tn u PrOW ' b "" f ° Und to be ■««* 

„„ S ‘ arlm? ” th , BUCk * n%9) analysis of h '? h - Pressure seals, followed by Allaire’s 11972) 
eccentric seal analysis and Childs’ ('19831 Hirs’ h.,Ur . , - y iaire s v iy7 ^) 

° r M " u,M z: * ■; 

Sharrer “V' ' ( C r dyn ‘‘ miC Coefflci ™‘s »as firs, considered by 

with wIvTp^tKL 1 .^ ~ •; ~ of a seal 

of a polynomial. Their analysis ronfirmAH a , , , a c ^ eara nce function in the form 

to a change in the seal profile Similar resultTfor this" 86 *” rotord - vnamic coefficients due 
nqqn „ oi „ w „ • , , similar results tor this case were reported by San Andres 

“ d Nd ’°" <«">■ —d a similar 

seal Detailed theTmolkstk ^tudieTw JeveSted dS?”” al ° ng ‘ r* ° f the 

direction and a similar distortion occurs alone the drc„mf ? not hmjted to axial 

distorted seal profile is shown in Fig 1, The ckar^ces f^hU^r m°' A " f ample of a 
a thermoelastic analysis. r t * us P ro ^ e were obtained from 

profi™’ rhrX'rarTseal pmfiTlfbaT’ ^ *•> with arbitrary 

Si 5 Thu ^ « « W-^hTpmpUX Nalso y „ SguveST 

asataibn„£^ 

for'co*mputing t resereffidenrirectiy°in'^glo^) r t^ tn ^^^^^^ofU^MaTand^n^hod 

addition, a new procedure for ^comS s.a^o.ffi ,T rd ‘? , ‘ e system is P res ™«d. In 

is also discussed computing seal coefficients based on external load specification 

discussed abOTe^t^ I ro”ed I mlalysis > Is a emDl* Pt ’d a t ; i,h MiaI «™*™. is 

of a SSME TWbopumpZlihe r«SsIre ro™ ^ ° ‘ dis,or,ed i«"»age seal 

inlet and exit clearances compared to those of a similar seal with average 


THEORY 


Bulk Flow Governing Equations 


Mass conservation and force equilibrium considerations in the axial 
directions for the control volumes in figures 2a and 2b vieM rh* f ll 
tmuity, axial momentum and circumf^ntiallmtn^'^"' 


circumferential 
bulk flow con- 
incompressible 


1 9(hv) d{hu) dh 

R diS ~ + Tt = 0 ( 1 ) 
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h dp 

JR33 


dv v dv dv 

+ m + “^ } 


+ 


f t - \/u 2 T I ’ 2 + /p- — ^ — - ^/u 2 + (v - U )) 2 


(3) 


where the friction factors /. and f T are defined for the Hirs and Moody friction factor models 
in the appendix* 

The boundary conditions at the inlet and exit of the seal are given as, 


Poi — Po(0. J) — (1 


0 ^ 4 ( 0 ,^) 


(4) 


i>o(0,/3) = P 3r x 

(5) 

Po{L.3) - POr. 

(6) 


where p 0< and po,. are the entrance and exit pressures respectively, * is the entrance loss 
coefficient and psr is the pre-swirl ratio. 


Film Thickness 


The expression for film thickness h(:,d) as a 
fixed coordinate system, instead of a minimum 


function of eccentricity is derived in a 
film thickness” coordinate system. The 
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r igure 


^lnerentiai Fluid Volumes, 


/ continuity b) Momentum 


coordinate system (*,*) shown in Fig. 3 is fixed at th. t 

,s onented P^allel to a user defined-global coordin t e «entric position and 

operation of a uniform profile seal such as a str^h arV™ , Typ , ica11 * for acentric 

coefficients are computed as a function of eccentricifv in a ape j Cd sea1, the rotordynamic 

the hne °f minimum film thickness. The use of these copffi CO( * d,nate astern aligned with 

a stability analysis requires their transformation intn tn fficient j 30 application such as 

This procedure, which is valid for a seal with uXId Sp' ^ COOrdmate system. 

non-uniform profile in the circumferential dLcT/on Such 1 “* applicable for seals with 

these coefficients m the user defined coordinate system dirlft] the imputation of 

with the angle of minimum film thickness fanirle nf • . aS tbese ^efficients vary 

The cpal tmcKness (angle of eccentricity). J 

C specifies £ I* *""" “ A 

seal profile will be non-uniform if c varies with 0 ^ " eal and 50 on ' The 
eccentricity, is derived as a function of c(» 3) and th ■ J<dtness ’ w ^ich varies with 

the film thickness and its gradients are .hei bel T*™* E ' The «Wk» for 
specifying the film thickness in a fixed coordinate system tW “ !° Fig ' 3 ' Besides 
accurate, particularly at high eccentricities than the mnr 8en * T ? ex P ression »* more 
form, h 0 = c - E x costf - E ;/ sin3. ’ * m ° re com monly used approximate 


= - (E.cosj + E y ,inJ) - g 

^ = (£± c >§ - - Eyco^^wi + E ,sinJ) 


(?) 


\J(R + c ) 2 - (E^inf 3 - E y cos/ 3) 2 + ( E * stn & ~ Eycosfl ) ( 8 ) 


dh o 

~d7 


(R + c) dc 


— 

sJTR + c ) 2 - (E x sin 0 - E, 


yCOSj ) 5 


( 9 ) 
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Figure 3: Diagram for Deriving General Clearance Expression 


Solution Procedure for Zeroth-Order Equations 

Th* solution for zeroth order equations involves the direct integration of the three cou- 

sSk 

A Zt hod bv which the coupled partial differential equations are reduced to coupled 
ordinary 5^^* " ith res P ect 

present Analysis, ^simpler method based on cubic sphne^ isns«h T ^ s ™^od ^more 

SSJSX i IS. Si the elliptical sea. in Figure 5 represents a special case of the 

“ b ^.lS^!^ti. F 4SiJn. are arranged in the following fashion and integrated 
from inlet to the exit. 

' ^ 1 [ g u ( u o,vo,Po,^0^^0<^) 

& = g v (u o,t-o,Po,^s^.^) (10) 

^ . ^(uot VotPo.Tjf-t^'^) . 

' T1, Vd»^ing n « eaA’cfmul^eremfS^c^ion^^th^direction of the corresponding 
paTient! ?or A, n««^e“ In other words, a, the ,-th ttxial step, the ctrcumferenttal 
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gradients are computed using the values of u . , 

the previous step, in the current anal • and Pou-i.j) •«■<?., the values at 

is used to compute these gradients. NelsonTndXuv^ 110 !! I?®™* b . ased on cubic s Pl»nes 
the above numerical integration For the current 8 / u . d th * simple Euler ’ s method for 
and 5th order Runge-Kutta method as well ^re^ 4th 

First Order Equations 

abouflheTt^df st 0 ate fc eccent d rk “h^oUowtn/ ° f the rotor 

P = P0 + rn , U = U Q + euj and „ = , 0 + ^ mg eXpressions - h = h 0 + <h u 

order terms yields the^onoXgS order^qTa^ns. 1 ' 3 negleCting second and higher 


dui h 0 Oi>i dh x 

°«r + * aj + «?“' + 


J_dh, 

R d/3 


dhi 

~dt 


u l - Uo 


dhi 

!h 


, & u \ . h 0 dp { 

°«r + 7ar + ' , °“° 


— _ / & U Q 1 <9t>o 

R dfS aT + R'df3^ hl 

dux , h 0 i'o <9ui 

<9- ~RT~d0 + ^ uUl + A '< v ' ~ A hhx 


/i 0 — + + /, „ V’odu, 

pi? dj3 + ft ° U °17 + + Buux + B v vx = B h h x 


(H) 

( 12 ) 

( 13 ) 


“**• wiabUs *. *. * 

both the Hir’s and Moody’s fricUon factors models "* gmn m the a PP endi * for 

The boundary conditions for the fir st order solution are (Ne,s„„ and Nguyen, 19M) , 
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Pi(0,id) = -(1 + Op u o(°- J ) ui(0,^) 

(14) 

r* 

o 

Xju 

II 

o 

(15) 

Px(LJ) = 0 

(16) 

Assuming that the rotor whirls about its equilibrium position in an elliptical orbit whose 
sen^ma^^and semi-minor axes are ,V„ and y„ respectively then the pos.t.on ot the center 
of the rotor relative to its static eccentric position is given by, 

X = XoCOSQ 

(17) 

Y = Y 0 sina 

(18) 

where a — ujt and u> is the whirl frequency. 

Let Ax - £a f and Xy = ^ where c 0 is the nominal clearance, and; 


(p X - Atj-Plx + &tyPly 

(19) 

fu\ — Ae^nix T AfyUir/ 

(20) 

€Vl = Ae r Vix + Ae y Uiy 

(21) 

ehl = Afx^lx + Afy/lly 

(22) 

hi x = -c 0 coaacosfi 

(23) 

hiy = -C 0 sinasinl3 

(24) 

Assume a solution of the form: 


p lx = a x {z,i3)cosa + a 2 (z,0)sina 

(25) 

u lr = a 3 {:,0)co3a + a A {z,i3)sina 

(26) 

V\ x = a 5 (c,/3)cosQ + U6( = 

(27) 

p Xy = b l {zJ)co3a + b 2 {z,l3)sina 

(28) 

u Xy = b 3 (z, i3)cosa + b A {z,(3)sina 

(29) 

v Xy = b s {z,i3)cosa + b 6 {zJ)sina 

(30) 


Using the above substitutions in the set of first order equations yields 12 coupled linear 
partial differential equations. The same solution procedure that is used for the zeroth 
solution is used to numerically solve for variables a, and 6.. 
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The first order boundary conditions 


are expressed in the assumed solution variables as: 

a^O, J) = -(1 +{)pa 3 (0,J) 


a 2 (0.^) = -(i .f ^)pa 4 (0.d) 

«s(0,^) = 0 

MM) = 0 

<*i(M) = o 


MM) = 0 

The original equations involving 6,’s. 

two auxiliary functions of the form ° b<? llarmoruc and separated them i 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 


into 


ji\~jcosp gi(z)$infj 


here fi and g t arc assumed not to vary with /? j xt 

a second separation of variables substitution to^the first ^ hT" 1988a) thereby apply 
N-16). While the above form of assumed so.ution ( ** 

experimental results, an examination of „ . . 1 that agree w,th available 

revealed a * dependence, ZuonZ v 1 ! T ^ ° f fu " Cti - s /<(-*) and ft(c) 
circumferential gradients should therLre irnDroTth* ( °' 5) ' The inclusion of the « 
The a, and 6, in the current analysis are totalW* S ° Utl ° n f at h, S her eccentricities, 
thereby avoids the mathematical contradiction d^/ functions of ; and 3 which 

cases the results of the current ^ • ^ rthermore - in many 

then the earlier results. PP better a g reei nent with experimental results 

The solution procedure for the 12 linear 'Dnt?> ■ ■ 

order solution. The solution is performed with V th and that ° f the zeroth 

and also with a predictor-corrector method * Brnh^rt. j h Run S e -Kutta method 
».th Runge-Kutta based method betng the fastest ‘* lmost iden,ical results ' 


Dynamic Coefficients 

position is given by int^at2”fhe 'oldlr pwssMe * SUtic CCCentric 


- &F r 


-af , ; 


~ f [ epiCOsdR d3 dz 
Jo Jo 


; - 1 r (P\sin/3R dfj dz 


(38) 


(39) 


cients. In this equation, IndY 'definT t he° mlaU ve'dfspl’ th f r . 0t0rdynamic coeffi - 

are the components of the force due to first order pressure Sd"“ ^ ^ F " Fy 
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( 40 ) 



The original analysis discretized the circumference into a number of strips and the 
function values (/;,*) are assumed to be independent of 0 over each strip. The 
method improves this approach by allowing the a t and 6; to vary over each strip m ob g 

thC Substftu^e'eqs. 0 28-33 into 22-27 and in turn substitute the results into eqs. 43-45. Also 
substitute eqs. 20 and 21 into 43-45. This yields: 


K xx - -V/xxuf 2 = — / / aiCosflR dP dz 
Co Jo Jo 

(41) 

c xy ^ = - [ L r bxcos&R d3 dz 
* Co Jo Jo 

(42) 

- fry* + rn yT u> 2 = - / 1 a\sinpR dp dz 

y Co Jo Jo 

(43) 

1 f L fir 

C vy ^J = — / / b\3inpR dp dz 

yv C 0 Jo Jo 

(44) 

- C xx *> = — / f a 2 cospR dp dz 
Co Jo Jo 

(45) 

L J o 2 ' 0 dz 

(46) 

c ux c *; — — / / a2sinf3R dj3 dz 

y c 0 Jo Jo 

(47) 

K « - m »“ ! = 7 0 Jo L h ‘ ia0R dt> ds 

(48) 


These 8 equations are evaluated for at least two whirl frequencies to obtain Rations 
for the 12 dynamic coefficients. A least squares approach is employed for this step. The 
2D Integration performed numerically are an improvement over the average value approach 
employed by the previous researchers. 


Dynamic Coefficients based on External Load Specification 

In some cases, it is possible to specify the angle at which external load is supported by 
thesealTuring the operation of the turbomachine. This external load ,» equal and opposite 
to the resultant seal force. A new method of computing the rotordynamtc coefficents based 

■*« is loc>,ed i,era,i,el ? such * at th ' re ; s e .r- 

librium between P the external specified load and the resultant seal force. The angle at which 
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the resultant seal force acts is forced to alim n«r>0\ L 

For example, unit 3-01, an experimental seal under iesign ItVIsA^ angle ‘ 

ater suppom thE ,o * d * * ~ of 2 9 o» 

Determination of Steady State Force Equilibrium P„ s i,i on 

position. Atth^teady ^“uilibJh^po.'S,” ' W ° < ' ime " sions ,0 locate the operating 

ft — F x - Wsinfi ~ 0 
fy — Fy - IVcos/J = 0 

The modified 2-D Newton-Raphson iteration procedure is described below 

- 


Ax 


Ax 


*u. 

dx 

K, 




a Ui.jfi + fx\r.,.y, = 0 


9fx 

dy 

±y d A 

V dy 


(49) 


dx T Qy + fy\x„y, — 0 (gg) 

The gradient ^, Z ^ZT ^ T ^ ° f r ° t0r P ° sition (**.»).- 

This iterative fs'r ,“ at,1 “ t/TeT *“7 ^ ^ <*•«>• 

resultant seal forces. Once this equilibnlim noilH^^ “ temal , load ■« balanced by the 
proceeds as before. P lon 1S attam ed, the remaining analysis 

Verification Case: Allaire, et. al. 

preach results to the ir"?* 1 N 8Uyen-Nelson ap. 

similar direct stiffness, damping and cross ^coupled A “ ' hree app ™^ stl ™ 

figures 6, 7 and 8 respectively. P d stlffness vs - eccentricity as seen in 


seal length 
rotor radius 


^ eal 


Co 

fluid 

density, p 
viscosity, p 
A P 

rotor speed 
friction factor 
relative 

roughness, e/2c 0 0.000001 (stator) 
pre-swirl ratio 0.1 

■inlet loss, £ 0.5 


40.6 mm (1.60 in) 

39.9 mm (1.57 in) 

0.14 mm (0.0055 in) 

0.14 mm (0.0055 in) 

0.14 mm (0.0055 in) 

L02 

57.657 kg/m 3 (3.60 lbm/ft 3 ) 

7.4396x10-® Pa-s (1.5538x 10" 7 lb-s/ft 2 ) 
7.26 MPa (1050 psi) 

23700 rpm 
Moody’s 
0.0 (rotor) 
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Difwct D*ir*P ,n fl< G** (KN/nO 



Figure 5: Direct Stiffness .vs. Eccentricity for Allaire Example 



Figure 6: Direct Damping .vs. Eccentricity for AUaire Example 
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EXAMPLE OF AN ARBITRARY PROFILE SEAL: AN ELLIPTICAL SEAL 

The above analysis for an arbitrary profile is applied to the case of an elliptical seal with 
axially varying curvature. The results for a similar linearly tapered elliptical seal were 
initially reported by San Andres (1992). The motivation for this study is two fold The first 
is to show the general steps involved in the analysis of arbitrary profile seals and the other 

15 TwoT^ T T t6rmS ’ effect of a ^change in profile on the dynamic coefficients. 
rh _ nth of curvature are considered for this analysis: one with a linear axial profile and 
the other a quadratic axial profile. For this study, the mid-point clearance of the quadratic 

profile is made 7 5% of (c; + c B )/2, i.e., 0.75 times the mid-point clearance of a Unear profile 
with similar inlet and exit clearances. 

The equation of an ellipse is given by, 


x — acosfi 
y = bsinfi 


(51) 

(52) 


where a and b are the semi-major and semi-minor axes respectively. At any angular 
position l3 along the circumference, the radius r of the eUipse is given by, 8 


r = \J{acosii ) 2 -f (6sinJ) 2 
and the clearance c at this location is given by, 

c = r - R 


(53) 


(54) 


where R is the radius of the rotor. 

thli h :; e 7;r^ 8e f ni - minor axes of the eUi P se * some functional form along 
the length of the seal, the clearance is given by, 6 


< z i0)=\/(fi{:)co3i3) 2 + (f 3 (z)aini3) 2 - R 


(55) 

where fi(-) and / 2 (-) are the semi-major and semi-minor axes variations along the 
z-axis. The gradients of this clearance function are given in the appendix. 

The ellipticity 6 is defined as (Fig. 5), 




(56) 


where c x and c y are clearances at semi-major and semi-minor axes respectively and, 

c x — Ci at inlet 


c x — c r at exit 

and from above, 

c y = C x (l ~ £) (57) 

When 6 = 0, the elUpse reduces to a circle and for 6 = 1, the seal contacts the rotor 
The appendix provides the functions /,(z) and f 2 (z) for a Unear profile and a quadratic 
profile as afunction of delta . The results shown are for a centered seal as a function of 
ellipticity. The dynamic coefficients are normaUzed with respect to the coefficients for the 
linear profile case at 6 = 0. The values used for this normaUzation are K xx = 44975 kN/m 
(256883 lb/in), C xx = 21.78 kN-s/m (124.4 lb-s/in) and k xy = 15821 kN/m (90364 lb/in). 
The seal parameters for this case are given below. 
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Figure 7: Diagram for Deriving Elliptical Seal Clearance Expression 


Seal Parameters for Elliptical Seal 

seal length 

16.66 mm (0.656 in) 

rotor radius 

48.39 mm (1.905 in) 

c t 

0.069 mm (0.00273 in) 


0.099 mm (0.00390 in) 

Co 

0.069 mm (0.00273 in) 

fluid 

L02 

density, p 

1041.7 kg/m 3 (65.03 lbm/ft 3 ) 

viscosity p 

129.6x 10 -6 Pa-s (0.188x10-® lb-s/ft 2 ) 

A P 

25.39 MPa (3681 psi) 

rotor speed 

22700 rpm 

friction Factor 

Moody’s 

relative 

0.0 (rotor) 

roughness, e/2c 0 

0.03 (stator) 

pre-swirl ratio 

0.2 

inlet loss, £ 

0.33 


The plot for direct stiffness (Fig. 9) shows the effect of a small change in profile on the 
direct stiffness. For the linear case, there is a complete loss of stiffness at around S = 0.65. 
The stiffness for the quadratic profile is almost twice that of the linear profile. Also, it 
retains its stiffness over a much wider range than the linear profile. The difference in the 
other coefficients (Figs. 10,11) are relatively small. 

CASE STUDY OF A DISTORTED SEAL 

The distorted clearance profile for an interstage seal of the Space Shuttle Main Engine 
High Pressure Oxygen Turbopump (SSME- ATD-HPOTP) is shown in Fig. 1. The distorted 
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Normalized Olraci Dampmg. Cxx, Cyy Normalized Direct StHtnaaa. Kxx. Kyy 



Figure 8: Normalized Direct Stiffness vs. Ellipticity (Elliptical Seal) 



Figure 9: Normalized Direct Damping .vs. Ellipticity (Elliptical Seal) 
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clearance profile of this seal is obtained from a thermoelastic analysis. The clearances are 
provided at six axial planes along the length of the seal with 68 clearances at each plane. 
The clearances along the circumference are located roughly equidistant. 

The rotordynamic coefficients of the distorted profile are compared with those computed 
using average clearances at inlet and outlet respectively. The geometry and operating 
conditions at full power level FPL are given in the following table. 


Seal PaiXLTneters for Distorted Seal Unit 3-01 

seal length 

16.66 mm (0.656 in) 

rotor radius 

48.39 mm (1.905 in) 

avg. c, 

0.149 mm (5.87 mils) 

avg. 

0.148 mm (5.81 mils) 


0.149 mm (5.87 mils) 

fluid 

L02 

density, p 

1041.7 kg/m 3 (65.03 lbm/ft 3 ) 

viscosity p 

I29.6x 10 -6 Pa-s (0.188xl0~® lb-s/ft 2 ) 

± P 

35.25 MPa (5112 psi) 

rotor speed 

25000 rpm 

friction factor 

Moody’s 

relative 

0.0 (rotor) 

roughness, e/2co 

0.8518 (stator) 

pre-swirl ratio 

0.2 

inlet loss 

0.3 


The distorted seal profile is fitted with bi-cubic splines. The purpose of this spline fitting 
is two fold; one is to interpolate clearances at any given axial and circumferential location 
and the other is to numerically compute axial and circumferential gradients of the seal 

profile at any required location. . , , , , , 

According to the manufacturer’s specifications, the side-load on the seal acts at a con- 
stant angle of 290°. The seal coefficients for this variable profile seal are computed as a 

function of side-load acting at this angle. . . . , 

Figure 12. shows the relation between seal forces and eccentricity. No load operation 
requires the seal to be slightly off-centered due to the distortion in the seal. Figs. 13,14 
and 15 show how the dynamic coefficients vary with externally applied load and the effects 
of distorted clearance profile versus average profile (average of clearances at the inlet and 
exit circumferences). The coefficients are seen to be sensitive to high loads and also show 
significant changes due to the distorted profile, i.e., see Fig 15. 

CONCLUSIONS 

The current approach has improved on the original Nelson-Nguyen method (NNM) by; 

(а) Employing a continuous interpolation of the first order variables in 
the circumferential direction, and 

(б) Utilizing cubic splines instead of Fourier series for the circumferential 
interpolation of both zeroeth and first order variables. 

In addition the current method models seals with arbitrary clearance profiles in the circum- 
ferential and axial directions. This capability was demonstrated with the operating seal 
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Direct Stiffness 


OB 


0.7 

0 6 

0 5 
o 

c 0.4 

s 

o 

LU 

03 
0 2 
Of 
0 0 

0 - <50 600 750 900 .050 ,200 ^ 

External Load (load angle 290} (lb) 

Figure 10: Ercen.rici.y vs. Preload for Opera.iug Seal Clearar.ee Profile (Fig ,.) 






Figure 12: Direct Damping vs. Preload for Operating Seal Clearance Profile (Fig. 1) 



Figure 13: Cross Coupled Stiffness.vs. Preload for Operating Seal Clearance Profile (Fig. 
1 ) 
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profile of a SS ME- HP OTP seal F' 11 

equilibrium po si .i„„ of the sea! gi™ ^ «- operating 
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APPENDIX 


The coefficient 

expressions for the first order equations are defined as 


A u = 

h 4 rr + + + (/.p.o + /.Pro) 

(58) 

A v = 

h 0 du 0 F, 0 , F rn 

R dp + UoV °ir ~ + u °( v ° ~ 

(59) 

Ah. = 

1 d Po du 0 v 0 dv 0 u 0 

p d z U °lh ~R ~d0 + ^ h "° Ul ' 0 + /lrf r °) 

(60) 

B u = 

h °ip + u o l ’07p- + u 0 (v 0 - w)^- 
U ~ 1 -0 U r 0 

(61) 

B„ = 

h 0 dvo 2 F,„ F ro 

RlB " Vo u7o + {v °- w) cr 0 + + frUr o 

(62) 

B h = 

1 <9p 0 5i'o t'o dv 0 tt h,o h 

pd0 U ° dz R~di3 + VoU '°-^ + ( v o-*)U Vo^ 

(63) 



(64) 


with further definitions for Moody's and Hirs’ friction models given in the following 


Moody's Model 

tftra' Model 

U T 0 = (uo 2 + (vg — U ;) 2 ) 1 ^ 2 
= (uo 2 + u 0 2 ) 1/2 
R r0 = ^&(u£ + (v 0 - w ) 2 )V 2 

= ^(u 0 2 + no 2 ) 1 / 2 
/rO= ^[1 +(10 4 ^ + lO 6 ^) 1 / 3 ] 
/ rt = o*w ll + (10 <^ +10 « H L ) ./3| 

9ro = 00 S 0 10 Vo 4 fe + ^)- 2 / 2 

^ = 00 s;° (10^ + ^)- 2 /3 

ftrO = T«( 10 4 f, + ^1)1/3 
= ”“ 5S (100^ + $■/■ 

/r — /rO 
/a = /»0 
/r ~ /ro/2 
/* — /«o/2 
5r — 9ro / 2 

5a = W2 

/i r = /i r o/2 
:= /i.fo/2 

jt 1 — A ~ 

r ro - 2 ~ 

JP _ /* -3# 

I 

UtO — (uo 2 + (uq — it) 2 ) 1 / 2 
^0 = («0 2 + Uo 2 ) 1 / 2 

RtO - ^-{ul + (l-o - w) 2 ) 1 ' 2 
R -0 = ^(V + i’0 2 ) 1/2 

frO = nr[^ff r0 j mr 
/,o = 

3 r 0 = -m r n r [fl r0 ] m ' 

9» o = 

h rQ - -m r n r [fl r0 ] m ' 
hjo = 

fr — frO 

ft = ft 0 

/r = /ro/2 
/* = /»o/2 
ffr = ffro/2 

= <7*o/2 

h r ~ ffro/2 

= W 2 
jr = 
x ro 2 
r — f»-9» 

r SO — 2 


The first-order governing equations are expressed in terms of the a, and 6, functions as 
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h Q da\ dho uodho hov^Oa^ hoiiQ da$ 

' + (^-«o^)a3 + W 4 (.4,,-- — )a 5 -- — + 


p dz 


,, . ,C?U 0 ,1 <9l’o s V , UqVq . 

c 0 [(>U + «o(^T + - —jf sin t\ 


( 65 ) 


ha da.2 , , . 5/io. , . uodho . , h 0 voda A 

-—- h 0 u>a 3 ~ (A u - u 0 -^-)a A + (A v - —-^-)a 6 = c 0 ^u 0 cosJ ^~~dJ 

/iqUo 5ag 


5/3 

5as dh 0 1 dh 0 r 5uo 1 5t’o, 3 l ’o . 21 da 5 

'“”87 + 87° 3 + R M as = Col( 77 + W M - fi”"" 1 “ R 88 

da$ dho 1 dho ho da$ 

ho ~d7 + 57 a4 + R w a% = “ Cou,w - -ff aJ 

, 5a 5 „ , „ , *0 dai h 0 v 0 da 3 

/IqUq — *T 4" B v &% 4* Cqa3^C05,^/ ^ 

5a6 nLD ^° 2 ^° l ’° ^ a ® 

/l ° u °'57 + S,ia “ ~ hou;as + 5,,a6 “ ~~Rp~dd Wdd 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


hadbi ,. dh 0 ., , . , , , uodho . h 0 v 0 db 3 

— |- (A u — u<j-t — )^3 + houib A + (A„ — — — )b$ — —Cou>uosin0 

p dz dz 


R 5/3 


R 5/3 


+ 


h 0 u 0 db$ 
~R~~dJ 


(71) 


h 0 db 3 dh 0 ,. , , u 0 dh 0 

~~q 7 ~ ho*>b 3 + (A tl - u 0 -^)b A 4- (A„ - 

r , ,5u 0 1 di’o., . , u 0 i'o ,, h 0 v 0 db A h 0 u 0 dbo 

= + “»< 87 + r 8j >)s '^ + ~r~"”^ ~ —Jj3 ~lg 

, db 3 dh 0 , 1 dh 0 , . . h Q db 5 

h ° 87 + 87‘” + R 88 ts = C ‘ VJam ' 3 - R 88 


, V.-, -.-u, 1 dho, ,1’0 , , ,duo 1 dvo. , ,, hodbo 

Ro— + —(>4 + g Jjt. = £01-^00,8 + (-jr + R 88 - R 88 


564 5/lo 

57 + 57 


db* ho db\ /jni'o 56s 

*o«o^ + S„6, + SA + fco-i. = -J-pjj - T M 

5 fr« „ , , . _ , - ho db2 h 0 v 0 56 6 

/iqUq ^ ^ 4* B u &4 — floats + — —coBh$mi3 

The clearance functions for the elliptical seal are given in the next table as; 


(72) 

(73) 

(74) 

(75) 

(76) 
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linear taper 

quadratic curve 

fi(z) = a\ a 2 z 

/t(-) = <*i + «2- + < 23 ^ 

/2(~) = ^1 + ^2- 

A( - ) = ^1 + i>2~ “1“ 

= a 2 

^ - a 2 + 2a 3 z 

-© 

II 

— b 2 -r 26 3 z 

a! = ii + Ci 

a x = + Cj 

a 2 = r( c « - c i) 

^2 — Trfcr. — ^Crn 4" 3Cj) 
■" p(Q 2 c t „ 4- Ci) 

&i = + ( 1 - *)c; 

b 1 = /? + ( 1 — ^)cj 

^2 = r(l - <5)(Cf: - Ci) 

^2 = -jf(l - ^)( c r - 4Cm + 3C{) 

*>3 = z( l ~ 6 )( C < ■ ~ 2c m + C»)) 


Gradients of the clearance function for elliptical seal are given by; 


c (--i) = (fi(:)cosi3) 2 + ( f 2 {z)sinf3 ) 2 - R 


(77) 


dc _ fif[cos 2 & + f 2 f' 2 sin 2 )3 

Tz' 

dc _ (/| - ff)cosl3sin/3 

^ \j{ficosj3) 2 + (f 2 sin0) 2 



(78) 

(79) 
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